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Abstract 



^ ■ We study the n — ► 77 , n — > 7/i /j + and 77 — > 7r°77 decays using an 

■ extended three-flavor Nambu-Jona-Lasinio model that includes the 't Hooft 

instanton induced interaction. We find that the 77-meson mass, the n — * 
77 , r? — > 7/x _ /i + and 77 — ► 7r°77 decay widths are in good agreement with 
the experimental values when the {7a(1) breaking is strong and the flavor 
SU(3) singlet-octet mixing angle 9 is about zero. The calculated 7777* 
Q H ' transition form factor has somewhat weaker dependence on the squared 

four-momentum of the virtual photon. The effects of the ?7a(1) anomaly 
on the scalar quark contents in the nucleon, the S^jv and T<kn terms and 
^ ■ the baryon number one and two systems are also studied. 

«?; 

1 Introduction 

It is well known that the QCD action has an approximate Ul(3) x Z7r(3) chiral 
symmetry and its sub-symmetry, Ua{1) symmetry, is explicitly broken by the 
anomaly. The Ua{1) symmetry breaking is manifested in the heavy mass of the 
7/ meson. The physics of the 77 and rf mesons have been extensively studied in the 
1/Nc expansion approach In the Nc — > 00 limit, the Ua(1) anomaly is turned 
off and then the 77 meson becomes degenerate with the pion and the 77' meson 
becomes a pure ss state with m^,(N c — > 00) ~ 2m 2 K — m\ ~ (687 MeV) 2 [0. So 
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the Ua{1) anomaly pushes up by about 400MeV and by about 300MeV. 
It means that not only the rj' meson but also the rj meson is largely affected by 
the Ua(X) anomaly. 

In order to understand the role of the £7a(1) anomaly in the low-energy QCD, 
it may be important to study the 77-meson decays as well as its mass and decay 
constant. Among the 77-meson decays, 77 — > 77 , 77 — > , ~fl~l + (I denotes e and/or 
fi) and 77 — > 7r°77 decays are interesting. They have no final state interactions and 
involve only neutral mesons so that the electromagnetic transitions are induced 
only by the internal (quark) structure of the mesons. 

The 77 — > 77 decay is related to the Adler-Bell-Jackiw (ABJ) triangle anomaly 
through the partial conservation of axialvector current (PCAC) hypothesis. 
One of the useful and widely used frameworks for studying the phenomena related 
to the axial-vector anomaly and the spontaneous chiral symmetry breaking is the 
chiral effective meson lagrangian given by Wess and Zumino 0] and developed 
by Witten |J. The 77, 77' — > 77 decays have been studied using the Wess-Zumino- 
Witten (WZW) lagrangian with the corrections at one-loop order in the chiral 
perturbation and it has been shown that the two-photon decay widths can be 
explained with the 77-77' mixing angle 9 ~ —20° ||. From the chiral perturbation 
theory (ChPT) M point of view, the WZW term is derived in the chiral limit 
and is of order p 4 . As discussed in ||, to reliably calculate 577(3) breaking 
effects of the 7/, 7/ — > 2 / y decays, the low-energy expansion to order p 6 has to 
be carried out. However in || full analysis of order p 6 has not been performed. 
Furthermore, because of the C/a(1) anomaly, the singlet channel decay amplitude 
Vo ~~ * 77 derived using PCAC + ABJ anomaly should be modified so as to become 
renormalisation group invariant ||. 

The 77 — > decay is closely related to the 77 — > 77 decay since it is 

considered as 77 — * 77* — > / yl + l~. By observing the muon pair invariant mass 
square spectrum of the 77 — > ^yl~l + decay, one is able to obtain the transition 
form factor for the 7777* vertex. It gives us the information of the size of the 77- 
meson. In |H| the transition form factors for the 7r°77*, 7777* and 7/77* vertices 
have been studied in the vector meson dominance model, the constituent quark 
loop model, the QCD-inspired interpolation model by Brodsky-Lepage and the 
ChPT. However none of the models have taken into account the effects of the 
Ua(1) anomaly explicitly. 
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Unlike the 77 — > 77 decay, there is no low-energy theorem for the 77 — > 
7r°77 decay and therefore it is not trivial. In ChPT |p"T|-|i3|l, there is no lowest 
order 0(p 2 ) contribution to the 77 — > 7r°77 process because the involved mesons 
are neutral. Likewise the next order 0(p 4 ) tree diagrams do not exist. Thus 
the 0(p 4 ) one- loop diagrams give the leading term in this process, but the con- 
tribution is two orders of magnitude smaller than the experimental value. This 
is because the pion loop violates the G-parity invariance and the kaon loop is 
also suppressed by the large kaon mass. The 0(p & ) contributions are dominant 
and the result is a factor two smaller than the experimental value. Although 
these results based on ChPT are not too far from the experimental value, it is 
noted that the higher order 0(p 6 ) terms in the perturbation expansion are larger 
than the leading 0(p 4 ) terms and the results contain ambiguous parameters that 
cannot be determined well from other processes. 

The purpose of this paper is to study the 77 — > 77 , 77 — > 7 l~l + and 77 — > 
7r°77 decays in the framework of a chiral quark model so that the quark structure 
of the 77 meson is explicitly taken into account. In such a model the explicit 
chiral symmetry breaking by the current quark masses can be included in a 
nonperturbative way. The effects of the {7a (1) anomaly can also be represented 
by the coupling of light quarks to the instanton as was pointed out by 't Hooft 
jnj]. We here take the three-flavor Nambu-Jona-Lasinio (NJL) model as the 
chiral quark model. The model involves the {7l(3) x Ur(3) symmetric four-quark 
interaction and the six-quark flavor-determinant interaction |1^ , [15| incorporating 
effects of the {7a (1) anomaly. It is widely used in recent years to study such topics 
as the quark condensates in vacuum, the spectrum of low-lying mesons, the flavor- 
mixing properties of the low-energy hadrons, etc. jn^ -j2jj. In this approach the 
explicit chiral symmetry breaking and the Ua(1) anomaly on the 77 — > 77 , 
?7 — > 'y Z~Z + and 77 — > 7r°77 decay amplitudes can be calculated consistently with 
those on the 77-meson mass, 77 decay constant and mixing angle within the model 
applicability. Furthermore one is able to study how the 77-meson properties change 
when the strength of the £/a(1) breaking interaction is changed in this approach. 

We have studied the 77 — > 77 decay p^] and the 77 — > 7r°77 decay [p3| 



111 



this approach and found that these decay widths are reproduced when the £7a(1) 
breaking interaction is much stronger than the previous studies in the three-flavor 
NJL model |T6|-pTf. The {7a(1) breaking six-quark flavor-determinant interaction 
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induced by the instanton [14]] gives rise to flavor mixing not only in the pseu- 
doscalar channel but also in the scalar channel. It is further argued that the 
instanton can play an important role in description of spin-spin forces, particu- 
larly for light baryons pI-2T)|. Since the Ua(\) breaking interaction is found to 
be rather strong, it is important to reexamine the effects of the Ua(1) breaking 
interaction on the scalar quark contents in the nucleon, the and T^kn terms 
and the baryon number one and two systems. 

The paper is organized as follows. In sect. 2 we explain methods for cal- 
culating the ?7-meson mass, mixing angle and decay constant in the three-flavor 
NJL model. We describe the calculations of the rj — > 77 , rj — > 7/i~/i + and 
rj — > 7T°77 decay amplitudes in sect. 3. The numerical results of the //-meson 
decays are presented in sect. 4. We study the uu, dd and ss contents in the 
nucleon and the and T^kn terms in sect. 5. Sect. 6 is devoted to the study 
of the effects of the {7a(1) breaking interaction on the baryon number one and 
two systems. Finally, summary and concluding remarks are given in sect. 7. 



2 ?7-meson in the three-flavor NJL model 

We work with the following NJL model lagrangian density: 



C 

Co 

U 



ip {id^ — m) ip . 
G, 



a=0 



G D { det [V>i(l - 75)^ + det ^(1 + 75)^ } • 



(1) 
(2) 

(3) 
(4) 



Here the quark field ip is a column vector in color, flavor and Dirac spaces and 
X a (a = . . . 8) is the U(3) generator in flavor space. The free Dirac lagrangian 
Co incorporates the current quark mass matrix rh = diag(m u , m<j, m s ) which 
breaks the chiral Ul(3) x Z7r(3) invariance explicitly. £4 is a QCD motivated 
four-fermion interaction, which is chiral Ul(3) x Ur(3) invariant. The 't Hooft 
determinant Cq represents the Ua(1) anomaly. It is a 3 x 3 determinant with 
respect to flavor with i, j = u, d, s. 

Quark condensates and constituent quark masses are self-consistently deter- 
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mined by the gap equations 

M u = m u - 2G S (uu) -2G D (dd) (ss) , 
M d = m d -2G s (dd)-2G D (ss)(uu), 
M s = m s -2G s (ss) -2G D (uu)(dd) , 

with 

rA d 4 P 



iaa) 



p^Y -M a + ie 



(5) 



(6) 



Here the covariant cutoff A is introduced to regularize the divergent integral and 
means trace in color and Dirac spaces. 
The pseudoscalar channel quark- ant iquark scattering amplitude 



(p 3 ,p 4 ;out|pi,p 2 ;in) = (2vr) 5 (p 3 + p 4 - p 1 - p 2 )% 



<N 



(7) 



is then calculated in the ladder approximation. We assume the isospin symmetry 
too. In the rj and rj' channel, the explicit expression is 



T - = 

■Lqq 



u(p 3 )\ 8 i'j 5 v(p A ) 

u(p 3 )^°il5v(P4,) 



T 



A(q 2 ) Btf) 
B(q 2 ) C(q 2 ) 



v(p 2 )\ 8 i-f 5 u(p 1 
v(p 2 )\°i'j 5 u(p 1 



with 



B(q 2 ) 
C(q 2 ) 



{2(G G 8 - G m G m )l\q 2 ) - G 8 } , 

2 

detD(g 2 ) {~ 2 ( G ° Gs ~ G m G m )T m (q 2 ) - G m } , 



detD 



— {2(G G 8 - G m G m )I 8 (q 2 ) - G } 



and 



Go = 
G 8 = 
G m = 



±G s -±(2(uu) + (ss))G D , 
±Gs-l((ss}-4(uu))G D , 



1 



3^2 



((ss) - (uu))G D . 



The quark-antiquark bubble integrals are 



I°(q 2 ) = 

/V) = 
I m (q 2 ) = if 



A d A p 
• A d 4 p 
A d 4 p 

WT 4 



X r ( c '/- D ) 



S F (p)\°ij 5 S F (p + g)A°ry 5 
S F (p)X 8 i^ 5 S F (p + g)A 8 ry 5 



TrM'V \S F (p)X°i l5 S F (p + q)\ 8 t l5 



(8) 

(9) 
(10) 

(11) 

(12) 
(13) 
(14) 

(15) 
(16) 
(17) 



with q = pi + p 2 . The 2x2 matrix D is 



nft ft_ Du(q 2 ) D 12 (q 2 ) 
{q) ~\D 21 ^) D 22 (q 2 ) 



with 



D^q 2 ) = 2G 8 I 8 (q 2 ) + 2G m I m (q 2 )-l, (19) 

D 12 (q 2 ) = 2G 8 I m (q 2 ) + 2G m I°(q 2 ), (20) 

D 21 (q 2 ) = 2G I m (q 2 )+2G m I 8 (q 2 ), (21) 

D 22 (q 2 ) = 2G I°(q 2 ) + 2G m I m (q 2 )-l. (22) 

From the pole position of the scattering amplitude Eq. ([8]), the ^-meson mass 
m v is determined. 

The scattering amplitude Eq. (|J) can be diagonalized by rotation in the flavor 
space 

/ u(p 3 )XHj 5 v(p 4 ) \ V lT / A(q 2 ) B{q 2 ) \ x 
qq ~ ~ \ u(p 3 )X°ij 5 v(p 4 ) ) le 16 { B(q 2 ) C(q 2 ) ) ±e 

Xle {v(p 2 )\^ 5 u(Pi) J ' { } 



u(p 3 )\^ l5 v(p A ) \ T ( D\q 2 ) 
u(jpz)X*'i l5 v(p A ) ) { D*(q 



2\ 



\ v(p 2 )X^ij 5 u( Pl ) ) ' [ZV 
with X 71 = cos9X 8 - sm6X°, X v> = sinflA 8 + cosflA and 

Ts=( C ° S6 ft " Sin / )■ (25) 
\ smv cosy I 

The rotation angle 9 is determined by 

tana) = cm^W) ■ (26) 

So 9 depends on q 2 . At q 2 = m 2 , 9 represents the mixing angle of the A 8 and 
A components in the r^-meson state. In the usual effective pseudoscalar meson 
lagrangian approaches, the rj and rf mesons are analyzed using the g 2 -independent 
rj-i]' mixing angle. Because of the g 2 -dependence, 9 cannot be interpreted as the 
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rj-i]' mixing angle. The origin of the g 2 -dependence is that the rj and rj' meson 
have the internal quark structures. 

For the 7/ meson, since the NJL model does not confine quarks, the 7/-meson 
state has the unphysical imaginary part which corresponds to the rj' — > qq decays. 
Therefore we do not apply our model to the 7/ meson in this article. 

We define the effective r/-quark coupling constant g v by introducing additional 
vertex lagrangian, 

= g^il^^^r, , (27) 

with X v = cos 8X 8 — sin 8X°. Here <fi is an auxiliary meson field introduced for con- 
venience and the effective 77-quark coupling constant is calculated from the residue 
of the gg-scattering amplitudes at the 77 pole, i.e. g> 2 = lim g 2_ tm 2(g 2 — m^)D v (q 2 ) 
The 77 decay constant f v is determined by calculating the quark-antiquark one- 
loop graph, 

fv = ^f 7 ^ Tr(C ' / ' D) ^7,75y^(p)7 75 A^ F ( P - q)} \ q , =m% . (28) 

One can easily show that in the C/a(1) limit, i.e., Gd = and m Uj d 7^ 777s, the 
77 meson becomes the ideal mixing state composed of u and d-quarks, namely, 
m,, = 777 w , g v = g n , f v = f n and tan 9 = -y/2. 

3 77-meson decays 

3.1 77 — ► 77 decay 

The P — > 77 (P = 7T°, 77, 77') decay amplitude is given by 

( 7 (A;i)7(A; 2 )|P(p)> = 1(2^ (h + k 2 - p) W e^^^ 2 CT T P ^ 77 (p 2 ) , (29) 

where e\ and €2 are the polarization vectors of the photon and the P — > 77 decay 
width r(P — > 77) is given by 

r(P-77) = |rp. 77 | ^. (30) 

In the analysis of the 77-77' mixing angle, the rj, rj' — > 77 decay widths and the 



current algebra formulae are used | 27[| . In order to see the assumptions used in the 
derivation of the current algebra formulae for the 77, 77' — > 77 decay amplitudes, 
we re-derive them here. 
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The starting points are the following PCAC relations which are modified by 
the ABJ anomaly. 

PAl = f^m^ + f^m^-^N^e^F^F^ , (31) 
3MJ = famftr, + fwml4v> ~ ^NJJ^e^F^F^ . (32) 



with 

(O|0M;|77) = faml , <0|3MJ|t/> = f^ml , 
(0\d»A°\ V ) = f 0v m 2 v , (0|c>M>') = /o^mj , 



(33) 



and £> a (a = 3,8,0) is defined as tr[{Q,Q}^] with Q = ±(A 3 + ^A 8 ). a is the 
fine structure constant of QED and F^ v is the electromagnetic field tensor. The 
simple 77-7/ state mixing is then assumed, i.e., 



It?) = cos#|r/ 8 ) -sinfllr/o) , ^ 
\rj') = sm9 \r]s) + cos 6* \rj ) , 



and 

(0\d»Al\ V8 ) = f 8 m 2 m , <0|<9M> > = / mJ o . (35) 

From Eqs. @-(|35|), f Sr) = f s cos6, f Sr) , = / 8 sin0, f 0r) = -/osinfl and f 0v < = 
fo cos 9 are obtained. The 77-meson field (f) n is expressed as follows. 

m;\ f s * fo V 

+ ^^J_ f ^ J_ _ ^ 
4vr M m2 ^ / 8 2^3 fo V 3 ^ 

Using the LSZ-reduction formula, the 77 — > 77 decay amplitude is given by 

T^ 77 = (j(h)"/(k 2 ) \rj(p)) 

= 1 J d A xe-^ x (U + m 2 v )( 1 (k 1 ) 1 (k 2 )\<j )v (x)\0} . (37) 

By inserting modified PCAC relation Eq. 

T"?7— i-77 — lim 

with 

,2 



r^M + rW,,^)] , (3 
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T (2) (2 
17— >77 \" 



x(7(/ci)7(fc 2 
(m 2 - p 2 ) 

x (i(hh(k2 



c °s# ,b sin6> nil - 
— PA} - —d»Al 

J8 JO 



(39) 



d xe 



-ip-x 



cos 9 1 



sin 9 



4 vr fc ^* T 



h 2VS fo 
0" 



(40) 



Using the current algebra techniques and from the consideration of the general 
structure of the matrix element, we can show that T« 77 (p 2 = 0) = | 
can be calculated easily, 



T(2) 

J -?7->77 



T (2) (2 
tj— >77 



z (2tt) 4 5 4 (fc! + A; 2 - p) e^^elKk 



a (rn* - p 2 ) 



a 1 / cos 9 ^sin6>\ 
VS V /a /o 7 

therefore the reduced invariant amplitude 7^_> 77 in the soft r\ limit is 
T 

-£ 71- 



'-ii—>*n(p 0) 



7T ^3 



1 ( COSd 9 /9 Sil1 ^^ 



Similarly, the reduced invariant amplitude % 



V-rr(P = 0) 



a 1 / 



sin# 



in the soft 77' limit is 
■cos#\ 



7T 



V fs 



2V2- 



h 

2 



(41) 

(42) 

(43) 
(P 2 = 



Finally, 7^ 77 (p 2 = 0) ~ 7^ 77 (p 2 = m 2 ) and 7^/_^ 77 (p 2 = 0) 
m 2 ,) are assumed in the usual analysis of the rj-rj' mixing angle |[27|| . Since 
and m^/ are not significantly smaller than the typical hadronic mass scale, the 
above results should not be taken quantitatively. Furthermore, because of the 
ABJ anomaly by the gluons (Z7a(1) anomaly) in the flavor singlet channel, the 
saturation of the flavor singlet axialvector current by the rj and rj field is rather 
questionable. 

In the NJL model, the on meson mass shell 7r°, 77 — > 77 decay amplitudes can 
be obtained by calculating the quark triangle diagrams shown in Fig. [I] and our 
results are 



T-k — >77 (jP 



mla) 



a 



71 

a 

~9v 



1 



cos^{5F; -2F'} 
-sin 0V2{5F% + F'} 



(44) 



(45) 
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l(h) 
7(h) 



Figure 1: The quark triangle diagram for 7r°, r] — > 77 decays. 



Here Fp (a = it, s and P = w , rj) is denned as 

f 1 , , 2(l-x)M a 
F a P = \ dx \ dy 



— mpx(l — x)(l — y) 

4M a . 2 / m P \ 

— arcsm — — — . (4o) 
mp \2Mj 



We can see that the integrand of Fp has an unphysical pole when mp > 2M a . It 
is due to lack of the confinement mechanism in the NJL model. 

In the chiral limit, the pion mass vanishes and F^ becomes 1/M U . In this 
limit, the Goldberger-Treiman (GT) relation at the quark level, M u = g^f^, 
holds in the NJL model and this leads to 7^o^ 77 = a/(irf n ) which is same as 
the tree-level results in the Wess-Zumino-Witten lagrangian approach It 
should be mentioned that we have to integrate out the triangle diagrams without 
introducing a cutoff A in order to get the above result though the cutoff is intro- 
duced in the gap equations in the NJL model. In the U(3)l x U(3)r version of 
the NJL model, the WZW term has been derived using the bosonization method 
with the heat-kernel expansion 1 2~9|, 30 1 . In their approach, 0(1/ A) term has been 



neglected and it is equivalent to taking the A — > 00 limit. 

3.2 Tj — > *y l~l + decays 

The P — > 7 l~l + (P = 7T°, 77 and / = e, fi) decay amplitude is given by 

(l(k)r( qi )l + (q 2 )\P(p)} = i(27r) 4 5 4 ( qi + q 2 + k-p) 

xe^ pa eu(q 1 )Yv(q 2 )^k p q a f P ^ l - l+ (q 2 ) , (47) 
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where q = qi + g 2 and u and v denotes lepton and antilepton's spinors. The 
lepton pair invariant mass square spectrum is 



dx 



a 3 (1 -xy ( 



96tt 



-m P 



x 



r 
x 



i 

2\ 2 



1 + 



2x 



Tp^i-i+^px) , with (1 > x > r 2 ) . (48) 

Here x = q 2 /m 2 p and r = 2mi/m P and the total P — > 7 / _ / + decay width is given 
by 



r(p^ 7 rr 



a 3 f j 

m P I ax- 

96tt 2 p Jr* 



[l-xf ( r^( | r 2 - 



Tp^i-i+impX) . (49) 



By calculating the diagram shown in Fig. we obtain the 7r°, 77 — > 7/ / 
decay amplitudes as follows. 

a: 



71 

a 

7i 



: 9^Gl(q 2 ) 



9 V 



3V3 



co S e{5G u v (q 2 )-2G s v (q 2 )} 
-smeV^^G^ + G^q 2 )} 



where 



G a P {q 2 ) = f 1 dx f 1 dy 
Jo Jo 



2(1 - x)M a 



4M„ 



(m P — q 2 ) 



Ml — mpx(l — x)(l — y) — q 2 x(l — x)y 

2 ( Q 



. 2 imp 
arcsm — — 
2M a 



arcsm 



\2M a 



(50) 



(51) 



(52) 
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with q = \[q*. We can find the following relations, G P (q 2 = 0) = Fp and 
1~P^,i-i+(q 2 = 0) = Tp^ 77 . 

From the observed data for the two-photon transition 77* — > P and the lepton 
pair invariant mass spectrum of the P — > decay, one can obtain the P77* 

transition form factor /p 77 * (g 2 ) defined by 

/P77* ) = ^ 7^ ^ • ( 53 ) 

T P ^i-i+{q l = 0) 

For the spacelike q 2 , G P (q 2 ) is given by 

with (3 = Jl — AM 2 /q 2 . We introduce the slope parameter Ap by 



A P " dq 2 ^ 11 '^ 



6 



(55) 



and Ap corresponds to the pole mass if one fits the (^-dependence of fp- ry *{q 2 ) 
by means of a single-pole term. In the case of the charge form factor / c (<? 2 ) of 
the charged pseudoscalar meson, 6/A 2 , = 6df c (q 2 )/dq 2 = (r 2 ) is the mean square 
charge radius of the charged pseudoscalar meson. Therefore, it is natural to 
consider that Ap is related to the size of the neutral pseudoscalar meson P. 

A n can be calculated easily in the chiral limit. Using Eqs. (0), Q52J ) and Q55j), 
we get A,,. = \/l2M u . On the other hand, the NJL model predicts A c = 2\/2f lv for 



the pion in the chiral limit |yj . Since A n is expressed in terms of the dynamical 
quantity of the model: M u in contrast with A c which is expressed in terms of the 
observed quantity: f^, the P77* transition form factor may be more sensitive to 
the dynamical structure of the pseudoscalar meson than the charge form factor. 

3.3 77 — > 7r°77 decay 

The i] — > 7r°77 decay amplitude is given by 

(^(pMh, e 1 ) 1 (k 2 , e 2 )\ V (p)) = 2(2vr) 4 5 4 (^ + k x + k 2 - vVx^u • (56) 

The dominant contributions to this process in this model are the quark-box di- 
agrams given in Fig. |3|. Following the evaluation of the quark-box diagrams 
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Figure 3: The quark box diagram for rj — ► it 77 decay. 



performed in [32], we obtain 



with 



T, 



u, 



Hi/ 



u, 



u, 



4 
5 



u: 



fJ,U 



—i—j=(cos6 — \/2 sin 9) e 2 g v g 7 , 
a/3 



= tr|7 5 

x 7„ 



1 



aq Yu l 



1 



- M + ie 
1 



75" 



J-fa-fo- M + ie 
1 ] 

:7m- 



1 1 



tr j 75 ^ - M + ze 7f V + ^2 - M + ie 
1 1 

X lv^—L U 77 , - 



= tr|7 5 

x 7m 



+ - #1 - M + ze + - M + ze 
1 1 



- M + ie $ + #2- M + ie 

1 1 
-75- 



Ul,(h ~ fc 2 ), 
[7^ <- fc 2 ), 



(57) 



(5? 



(59) 



(60) 

(61) 
(62) 
(63) 



Here M is the constituent u,d-quark mass. Because the loop integration in (|57|) 
is not divergent, we again do not use the UV cutoff. Then the gauge invariance 
is preserved. The inclusion of the cutoff that is consistent with the gap equation 
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will break the gauge invariance and make the present calculation too complicated. 
Note that the strange quark does not contribute to the loop. 

On the other hand the amplitude T Mi , has a general form required by the gauge 
33| 



invariance 



A(x 1 ,x 2 )(k'[k$ - h ■ k 2 g 



+ B(x 1 ,x 2 ) 



-m 2 r} x 1 x 2 g^ v 



+ x x k^p v + x^kl 

mi 



(64) 



with X{ = p ■ ki/m 2 . With A and B, the differential decay rate with respect to 
the energies of the two photons is given by 



d 2 T 

dx\dx 2 



ml 



256vr 2 





2 " 


A+-B 




2 





mz 



2( Xl +x 2 ) + ^-l 



mi 



4x\X 2 



mz 



2( Xl +x 2 ) + ^-l 



mi 



1 2 



(65) 



Though the mass of r\ as a qq bound state depends on Gd, we use the experimental 
value m v = 547 MeV in evaluating (|65|) . The Dalitz boundary is given by two 
conditions: 
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(66) 
(67) 



In evaluating fl5"8|)-(|63|), one only has to identify the coefficients of p tJ- p u and g^ v . 
Details of the calculation are given in [^|. Defining A and B by 



+ 



we find A and B as 
A ■ 
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1 2 
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(69) 
(70) 



with 



(*i + fc 2 ) 2 
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2(x 1 + x 2 ) + 



mz 



- 1. 



(71) 



v v 
We evaluate .4 and B numerically and further integrate ( p3| ) to obtain the r] 

7r°77 decay rate. 
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4 Numerical results 



4.1 77- meson mass, mixing angle and decay constant 

We discuss our numerical results of the 77-meson mass m^, mixing angle 9 and 
the 77 decay constant f v in this subsection. The parameters of the NJL model are 
the current quark masses m u = rrid, m s , the four-quark coupling constant Gs, the 
{/a(1) breaking six-quark determinant coupling constant Gd and the covariant 
cutoff A. We take Gd as a free parameter and study 77 meson properties as 
functions of Go- We use the light current quark masses m u = = 8.0 MeV 
to reproduce M u = M d ~ 330 MeV (~ l/3Mjv) which is the value usually 
used in the nonrelativistic quark model. Other parameters, m s , Gs and A, are 
determined so as to reproduce the isospin averaged observed masses, = 138.0 
MeV, m K = 495.7 MeV and the pion decay constant f n = 92.4 MeV. 

We obtain m s = 193 MeV, A = 783 MeV, M u4 = 325 MeV and g nqq = 3.44, 
which are almost independent of Gd- The ratio of the current s-quark mass to 
the current u,d-quark mass is m s /m u = 24.1, which agrees well with m s /m = 
24.4 ± 1.5 (m = \{m u derived from ChPT |34|. The kaon decay constant 

fx is the prediction and is almost independent of Gd- We have obtained fx = 97 
MeV which is about 14% smaller than the observed value. We consider this is the 
typical predictive power of the NJL model in the strangeness sector. The quark 
condensates are also independent of Gd and our results are (wit) 5 = —216 MeV 
and (ss)5 = —226 MeV, which give 

_ ! m fK m K 



, , 1.01, — T1 JK — = 1.06. (72) 

—2m u (uu) —±{m u + m s ){uu + ss) 

As reported in the previous studies in the NJL model, the Gell-Mann-Oakes- 
Renner relations hold well for both the pion and kaon sectors. 

We define dimensionless parameters G e § = — Gd(A/27t) 4 AA^ ( ? and G s s = 
G r 5'(A/27r) 2 A' r c . The calculated results of the 77-meson mass m v and the mixing 
angle 6 are shown in Fig. [| and Fig. |^, respectively. When G e § is zero, our 
lagrangian does not cause the flavor mixing and therefore the ideal mixing is 
achieved. The "77" is purely uu + dd which corresponds to 9 = —54.7° and is 
degenerate to the pion in this limit. 

It may be useful to compare our results with those in the 1/Nc expansion 
approach. In the 1/Nc expansion approach, with the inclusion of the 0(1/Nc) 



15 



m n [MeV] 
600 



400 



200 




1.6 



Figure 4: Dependence of the rj meson mass on the dimension-less coupling con- 
stant G e jj. The horizontal dashed line indicates the experimental value. 



contribution by the Ua(X) anomaly in the flavor singlet-singlet channel, the square 
mass matrix of the mass term of the low-energy effective lagrangian in the rj^-rjo 
channel becomes as follows [01. 



4^2 



1™2 



M, 



3% - 3 m n 
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K 



(73) 



6 



YM 



(74) 



where 

x 2 

Here the isospin symmetry is assumed and YM means the pure Yang-Mills theory. 
The ghost field is defined by 

,2 



9 



(75) 



and G a is the gluon field strength tensor. By diagonizing the matrix given in 
Eq. ([73[) , we obtain 



m 



2 X 
m-ff H 



c , 



2m 2 K — 2m\ 



3N> 



c , 



O Y 

ok' (76) 
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Figure 5: Dependence of the mixing angle 9 on the dimension- less coupling con- 
stant G e jj. The horizontal dashed line indicates the experimental value. 



and 



tan# 



4 2 



l 2 2 
3< - m V 



(77) 



From Eqs. ( ff^) and ([fj]), it is obvious that in the Ua(1) limit the ry and rj' become 
the ideal mixing state with m v = m„ 



138 MeV and m v > 



2m\ 



mi 



687 



MeV. We compare the dependence of the r\ meson mass on the mixing angle 
calculated in the NJL model with that given in Eq. ( f7?|) in the 1/Nc expansion 
approach in Fig. ^| and find that the r] meson mass calculated in the NJL model is 
somewhat smaller than that in the 1 / iV^ expansion approach at the same mixing 
angle (except for the ideal mixing point) though the shapes are similar. The value 
of the mixing angle is one of the important quantity to understand the physics of 
the r] and r/ mesons. The mixing angle determined from the rj-r]' mass formula 
is often discussed in the literature. However since the r) meson mass is not so 
sensitive to the mixing angle as shown in Fig. it is not suitable to determine 
the mixing angle from the rj meson mass. 

We next discuss the rj decay constant f v . The calculated rj decay constant 
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Figure 6: Dependence of the rj meson mass on the mixing angle. The solid line 
indicates the result calculated in the NJL model and the short-dashed line in- 
dicates that in the 1/Nq expansion approach. The horizontal long-dashed line 
shows the experimental rj meson mass. 



is shown in Fig. [j]. It is almost independent of Gd and ~ f n . For example, 
f v = 91.2 MeV at G*§ = 0.7. Therefore it seems that the rj meson does not 
lose the Nambu-Goldstone boson nature though its mass and mixing angle are 
strongly affected by the Ua(1) breaking interaction. 

The charged pion and kaon decay constants can be directly obtained by mea- 
suring the 7r — > fiu^ and K — > fiv^ decays. On the other hand the decay constants 
for the light neutral pseudoscalar mesons tt°, rj and r\' cannot be obtained from the 
direct measurements. As discusses in Sect. 3.1, they are related to the P — > 77 
decay amplitudes in the soft limit. It is widely accepted that the soft meson 
limit is close to the on-shell point in the pion case. However it is the matter of 
discussions how good the soft r\ and 7/ limits are. We will discuss this point in 
the next section by comparing the rj — ► 77 decay amplitude which is directly 
calculated in the NJL model with that obtained using PCAC + ABJ anomaly in 
the soft f] limit. 
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Figure 7: Dependence of the r\ decay constant on the dimension-less coupling 
constant G e n ■ 



Since the NJL model does not confine quarks, the meson is not a bound 
state. Therefore we do not apply our model to the 7/ meson. On the contrary, 
the f) meson appears as a bound state. Nevertheless, one may ask whether the 
present model is applicable to the 77 meson since the binding energy of the r\ 
meson (2M U — m v ) is much smaller than those of the pion and kaon. In order 
to confirm it, we study the constituent u, d-quark mass dependence of the 77- 
meson properties. By changing the current u, d-quark mass m u d from 7.5 MeV 
to 8.5 Mev, the constituent u, d-quark mass M u ^ is changed from about 300 MeV 
to 360 MeV. Other parameters of the model are chosen so as to reproduce the 
experimental values of m n , m^, m v and /„.. This change causes the change of 
the binding energy of the 77 meson from about 50 MeV to 170 MeV. The change 
of the calculated i] decay constant is within 2% and the change of the calculated 
mixing angle is within 12%. This stability of the 77-meson properties indicates 
that the NJL model can describe the essential feature of the rj meson. 
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4.2 rj —> 77 decay 



The recent experimental results of the tt°,t] — > 77 decay widths are r(7r° — > 
77) = 7.7 ± 0.6 eV and T{r] -> 77) = 0.510 ± 0.026 keV || and the reduced 
amplitudes are 



f 

■ L r]— >77 



(2.5 ±0.1) x 10" 11 [eV]" 1 , 
(2.5 ±0.06) x 10~ n [eV]" 1 



(78) 
(79) 



Here we have used the two photon measurement result for the 77 — > 77 decay 
width. From Eq. fl4|) and Eq. (fH), we get T^ 77 = (5/3)7;o^ 77 in the £^(1) 
limit. Therefore in order to reproduce the experimental value of 7^ 77 , the effect 
of the {7a (1) anomaly should reduce 7^ 77 by a factor 3/5. 



►77 



2.50 x 



We first discuss the 7r° — > 77 decay. The calculated result is 7^o_ 
10~ n (l/eV) which agrees well with the observed value given in Eq. (|78|). The 
current algebra result is T n o^ J7 = a/(7rf n ) = 2.514 x 10~ n (l/eV), and thus 
the soft pion limit is a good approximation for tt° — > 77 decay. The chiral 
symmetry breaking affects 7^o^ 77 in two ways. One is the deviation from the 
G-T relation and another is the matrix element of the triangle diagram F(u, 7r°). 
Our numerical results are g n = 3.44, M u j = 3.52 and F(u,n°)M u = 1.015, 
therefore the deviations from the soft pion limit are very small both in the G-T 
relation and the matrix element of the triangle diagram. 

Let us now turn to the discussion of the 77 — > 77 decay. The calculated results 
of the 77 — > 77 decay amplitude 7^ 77 is given in Fig. |[ The experimental 
value of the 77 — > 77 decay amplitude is reproduced at about Crff = 0.7. The 
calculated 77-meson mass at G|f = 0.7 is m,, = 510 MeV which is 7% smaller 
than the observed mass. As for the effects of the symmetry breaking on 7^_^ 77 , 
our results are F(u,rj)M u = 1.36 and F(u,r])/F(s,r]) = 1.99 at G§ = 0.7. 
Therefore it seems that the soft 77 limit is not close to the real world in this case. 
Recently, Bernard et al. |36] calculated the 77 — > 77 decay width using a similar 
model. They used a rather weak instanton induced interaction and their result of 
r(?7 — > 77) is about 50% bigger than the experimental value. It is understandable 
from our analysis. 

G e § = 0.7 corresponds to Gd{ss) /Gs = 0.44, suggesting that the contribution 
from Cq to the dynamical mass of the up and down quarks is 44% of that from £4. 
In the previous study of the 77 and 77' mesons in the three-flavor NJL model JTB]-pT[], 
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Figure 8: Dependence of the rj — > 77 decay amplitude on the dimension-less cou- 
pling constant G e p. The horizontal dashed line indicates the experimental value. 



the strength of the instanton induced interaction has been determined so as to 
reproduce the observed rf mass though the rf state has the unphysical decay 
mode of the rf — > uu, dd. The strength determined from the rf is much smaller 
than G e § = 0.7, about 1/4 of the present case. One of the shortcomings of the 
NJL model is the lack of the confinement mechanism. It is expected that the 
confinement gives rise to the attractive force between quark and antiquark in the 
rf meson to prevent the rf meson from decaying to the quark and antiquark pair. 

In the PCAC + ABJ anomaly approach, if one assumes the SU(3) symmetry, 
i.e., fs = fn, and using Eqs. ((42]) and (f|3|), the mixing angle 9 and the meson 
decay constant in the flavor singlet channel fo can be determined so as to repro- 
duce the observed T(r] — * 77) and T(rj' — > 77). The results are 9 = —17.4° and 
fo/fir = 11. The pion and kaon loop corrections are then included in ChPT 
and the results are 9 = -21.8°, f 8 /f n = 1.25 and f /f n = 1.07. 

In the NJL model, the mixing angle at G c D ff = 0.7 is 9 = -1.3° and that 
indicates a strong OZI violation and a large (u,d)-s mixing. This disagrees with 
the "standard" value 9 ~ —20° obtained in the PCAC + ABJ anomaly approach 
and ChPT. One of the origin of the difference is that the mixing angle in the NJL 
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Figure 9: The flavor singlet decay constant f which is determined so as to 
reproduce the observed r\ — > 77 decay width as a function of 9 using the formula 
derived in the PCAC + ABJ anomaly approach. The solid line indicates in the 
case of fs = f n and the dashed line indicates in the case of f% = 1.25 

model depends on p 2 of the qq state and thus reflects the internal structure of the 
r\ meson. On the contrary the analyses in the PCAC + ABJ anomaly approach 
and ChPT assume an energy-independent mixing angle, i.e., 9(p 2 = m 2 ) = 6(p 2 = 
m 2 ,). Another point is that the £77(3) breaking of the rj — > 77 decay amplitude 
is rather large in the NJL model as we have shown above. It indicates that the 
soft rj and rj' limits used in the PCAC + ABJ anomaly approach are not so good. 
In the ChPT point of view, it suggests that the tree diagram contributions from 
0(p 6 ) terms may be rather large. 

If the PCAC + ABJ anomaly approach is considered not to be applied for 
the rf — > 77 decay, then one cannot determine the mixing angle 9 and the flavor 
singlet decay constant /o from the observed rj, rf — > 77 decay widths. However the 
relation between 9 and fo can be obtained from Eq. ( |42]) only by the rj — > 77 decay 
amplitudes assuming that fs is given. Fig. || shows the fo — 9 relations for the 
SU(3) value f 8 = f n and the ChPT estimation / 8 = 1.25 f n . One sees that a 
smaller flavor singlet component of the rj meson corresponds to a smaller /q. It 
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Figure 10: Dependence of the 77 — > 7/^~/i + decay width on the dimension-less 
coupling constant Crff. The horizontal solid line indicates the experimental value 
and the dashed lines indicate its error widths. 

is not strange that /o is smaller than f n since the rj' meson has larger gluonic 
component than pion because of the Ua(1) anomaly. 

4.3 fj — ► 7 l~l + decay 

The experimental value of the rj — > / ~ffi~fi + and rj — > r ye~e + decay widths 



and the calculated r] — > 7 / u _ / u + and 77 — > 7e~e + decay widths are shown in Fig. 

[To| and Fig. [Tl| respectively. At G<§ = 0.70 where the rj —>■ 77 decay width is 

reproduced, we obtain r (77 — > 7/i _ /j + ) = 0.407 eV, which is in good agreement 

with the experimental result shown in (0). On the other hand, our calculated 

4 We have used the two photon measurement result of the 77 — ► 77 decay width: T(rj — > 
77) = 0.510 ± 0.026keV. 



are |35j Q 




(80) 
(81) 
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Figure 11: Dependence of the r] — ► 7e~e + decay width on the dimension- less 
coupling constant Gff . The horizontal solid line indicates the experimental value 
and the dashed lines indicate its error widths. 



result of the r] — > 7e~e + at Gff = 0.70 is r (77 — > 7e~e + ) = 9.72 eV, which is 
about 50% bigger than the observed value given in fl5I|) . 

It is clear from the electron-positron pair invariant mass square spectrum of 
the 7] — > 7e~e + decay given in Eq. fl48|) that the strength is mostly saturated in 
the small region just above q 2 = (2m e ) 2 . Since in this region, the 7777* transition 
form factor f mi *{(f) is almost unity and the radiative corrections are found to 



be negligible [ftTH , the 77 — > ^e~e + decay width is strongly related to the 77 — > 
77 decay width and it is rather hard to explain the present experimental values 
of r(?7 — > 77) and T(r] — > r ye~e + ) simultaneously. 

The calculated 7777* transition form factor f m -y*(q 2 ) at G c ^ = 0.7 is shown 
in Fig. [12] and the calculated slope parameter defined in Eq. (fffi ) is = 0.94 
GeV. The recent experimental results of the slope parameter are as follows. The 



TPC/27 Collaboration at the SLAC gives A^ = 0.70 ± 0.08 GeV |8| and the 
CELLO Collaboration at the DESY gives A^ = 0.84±0.06 GeV So our result 
is somewhat larger than the experimental results. As for the slope parameter of 
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the 7r°77* form factor, our model gives the simple result in the chiral limit, 
i.e., A^ = \f\2M u ~ 1.126 GeV, which is also larger than the CELLO result 
= 0.75 ±0.03 GeV ||]. 

Ametller et al. |l(| studied the transition form factors for the P77* vertices 
with P = 7T°, 77 and 7/ using the most successful and/or traditional models of 
the low-energy QCD including the vector meson dominance model (VMD), the 
constituent quark loop model (QL), the QCD-inspired interpolation model by 
Brodsky-Lepage (BL) and the ChPT. They concluded that all the models con- 
sidered agree in the correct value for a mean Ap, but differ in the breaking pattern 
of the SU(3) flavor symmetry. Our approach is close to QL approach. In Ref. [[L(| 



a rather small constituent quark mass (M u ^ — 0.23 GeV) has been used and that 
is the main reason why QL prediction of A^ is in reasonable agreement with the 
experimental results. However, it should be noted here that since the constituent 
quark mass used in QL is smaller than m v /2, the unphysical imaginary part 
appears in the 77 — > 77* transition amplitude corresponding to the unphysical 
channel 77 — > uu, dd. 
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Our interpretation of the present result in the NJL model is as follows. Since 
the pseudoscalar mesons in the NJL model have the quark- ant iquark structures 
and therefore have the size. For the 7r°77* vertices, it is r w o ~ l/^/2M u ~ 0.43 
fm, which is about 67% of the experimental value r n o = ~ 0.64 fm. 

Our result of the ??77* vertex size is r v = \^6/A v ~ 0.51 fm, which should be 



compared with the experimental results r v ~ 0.69 fm |38j and r v ~ 0.58 fm [p9|| . 
The introduction of the quark-antiquark correlations in the vector channel may 
solve the problem of the difference of the sizes for the P77* vertices between our 
results and the experimental ones. 

4.4 j] — > 7r°77 decay 

The experimental value of the r] — > 7r°77 decay width is [|35| (] 



r (77 -> 7T°77) = 0.93 ±0.19 eV, (82) 

and the calculated 77 — > 7r°77 decay width is shown in Fig. At G c ^ = 0.70 
where the r\ — > 77 decay width is reproduced, we obtain T(r] — > 7r°77) = 0.92eV, 
which is in good agreement with the experimental data shown in fl82|) . 

In ChPT [11, 12 1 , since the lowest order 0(p 2 ) term and the next order 0(p A ) 
tree diagrams do not contribute to the 7] —>■ 7r°77 process, the 0(p A ) one-loop 
diagrams are the leading contributions. The calculated pion and kaon one-loop 
contributions to the r] — * 7r°77 decay width are T^(r] — > 7r°77) = 0.84 x 10~ 3 eV, 
t k\v -> tt°77) = 2.45 x 10~ 3 eV and T^ K {r} -> ^77) = 0.84 x 10~ 3 eV, which 
are more than two orders of magnitude below the observed width (^). This is 
because the pion loop violates the G-parity invariance and the kaon loop is also 
suppressed by the large kaon mass. At 0(p 6 ), there exists contribution coming 
from tree diagrams, one-loops and two-loops. The loop contributions are smaller 
than those from the order 0(p A ). In |[L1]| , the coupling strengths of the 0(p 6 ) 
tree diagrams are determined assuming saturation by the scalar meson ao(980) 
and the tensor meson 02 (1310) resonances as well as the p and u vector meson 
resonances. The obtained resonance contributions are r^ a; (^ — > 7r°77) = 0.18eV 
and Tp+ a;+ao+a2 (?7 — > 7r°77) = 0.18 ± 0.02eV. The contributions of other mesons, 



such as 6i(1235), /ii(1170), ^ii(1380) [0 an d other tree diagrams [jl3| , are found 



5 We have used the two photon measurement result of the r\ — > 77 decay width: Y(r\ 
77) = 0.510 ± 0.026keV. 
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Figure 13: Dependence of the 77 — ► 7r°77 decay width on the dimension- less 
coupling constant G e §. The horizontal solid line indicates the experimental value 
and the dashed lines indicate its error widths. 



to be small. Therefore, up to 0(p 6 ), ChPT estimation of the 77 — > 7r°77 decay 
width is about 1/4 of the experimental result. 

On the other hand in ||T| the 0(p 6 ) tree diagrams are evaluated by using the 
extended NJL (ENJL) model [^]. They calculated three contributions in ENJL, 
namely, the vector and scalar resonance exchange and the quark-loop contribu- 
tions. Their result is T(i] — > 7r°77) ~ 0.5 eV. They further introduced the 0(p 8 ) 
chiral corrections as well as the axialvector and tensor meson exchange contribu- 
tions, and finally obtained T(i] — > 7r°77) = 0.58±0.3 eV. A recent similar analysis 



in ENJL concludes somewhat smaller value for the decay width [43]. The dif- 
ference between our approach and that in |^1J are as follows. The ENJL model 
lagrangian has not only the scalar-pseudoscalar four quark interactions but also 
the vector-axialvector four quark interactions. However, the Ua{^) breaking is 
not explicitly included in their model and therefore the 77 — 77' mixing is intro- 
duced by hand with the mixing angle 9 = —20°. We stress that the introduction 
of the C/a(1) breaking interaction is important to understand the structure of 
the 77 meson. There is another difference. The coupling constants of the chiral 
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effective meson lagrangian predicted in the ENJL model are parameters of the 
Green function evaluated at zero momenta. On the other hand we evaluate the 
quantities at the pole position of the mesons. 

Calculated spectrum of the photon invariant mass square m^ 7 for the r\ — > 
7r°77 decay is shown in Fig [TJ|. As this spectrum is compared with those calcu- 



lated by ChPT in jHJ, we find ours to be similar to the one for d 3 = 4.5 x 10 
GeV -2 in |Tj| which involves an additional 0(p 6 ) contribution to the original 
lagrangian. Spectrum of the photon energy for the 77 — > 7r°77 decay is shown 
in Fig |l^, and given in ]lTJ in ChPT. Both are also similar, though there is no 
experimental result. 

In our calculation of the rj — > 7r°77 decay, we evaluate only the quark-box 
diagram given in Fig [| Since the vector and axialvector four-quark interactions 
are not included in our model, the only other contribution to this process is the 
scalar resonance exchange. In the ENJL model the contribution of the scalar 
resonance exchange is small JTJ]. The relevant process is 77 — ^ ao7r° — ► 777?°. 
Since the ao — > 77 decay amplitude is known to be small experimentally, the 
scalar resonance contribution is suppressed. Although the effect of the scalar 
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channel quark-antiquark correlation is not taken into account in our calculation, 
the scalar resonance contribution is (partly) included in the box diagram. 

If one includes the vector and axialvector four-quark interaction in the NJL 
model, the pseudoscalar meson properties are affected through the pseudoscalar- 
axialvector channel mixing and the model parameters with and without the vector 
and axialvector four-quark interaction are different. We expect that the models 
with and without the vector-axialvector interaction predict similar results for the 
processes involving only the pseudoscalar mesons with energies much below the 
vector meson masses. It is further argued that the contribution of the quark-box 
diagram to the 77 — > 7r°7r° process, that is similar to 77 — > 7r°77 , is quite close to 
that of the vector meson exchange in the vector dominance model [[44 |. 



5 Scalar Quark Contents in Nucleon 

The Ua(1) breaking six-quark flavor-determinant interaction £ 6 given in Eq. 
(U) gives rise to flavor mixing not only in the pseudoscalar channel but also in 
the scalar channel. Therefore, it is important to study the effects of the ?7a(1) 



29 



breaking interaction in the scalar qq channels. 

In the NJL model, it is known that the masses of the scalar mesons become 
just the twice of the constituent quark mass in the chiral limit. Introducing the 
explicit breaking of the chiral symmetry by the current quark masses pushes up 
the scalar meson masses above the quark-antiquark threshold. Since the NJL 
model does not confine quarks, we do not apply it to the scalar mesons in this 
article. 

In this section we will discuss the scalar quark contents in nucleon as well as 
the pion-nucleon sigma term E^-at and the kaon-nucleon sigma term T*kn- We use 
a rather naive additive quark model for the nucleon state, namely, the nucleon is 
made up by three noninteracting constituent quarks whose masses are determined 
by the gap equation shown in Eq. @. The scalar quark content of flavor a in 
the proton is then obtained as follows. 

(P\q a q a \P) = 2{U\q a q a \U) + (D\q a q a \D) , (83) 

with \U) (\D)) is the constituent u- (d-) quark state. The amount of the qq 
content of flavor a in a constituent quark Q of flavor b is deduced by using the 
Feyman-Hellman theorem, i.e., 

{Q b W\Q b ) = ^ (84) 

where m a and are the current and constituent quark masses, respectively. 
Here the isospin symmetry is assumed and therefore the following relations hold. 
(U\uu\U) = {D\dd\D}, (U\dd\U) = (D\uu\D) and (U\ss\U) = (D\ss\D). 

We have calculated the scalar quark contents in the constituent u-quark as 
functions of G e p and the results are shown in Fig. 16]. At G c p = 0, both (U\dd\U) 



and (U\ss\U) vanish as expected. Because of the larger mass of the strange quark, 
(U\ss\U) is strongly suppressed compared to (U\dd\U). At = 0.7, our results 
are 

(P\uu\P) = 3.30, (P\dd\P) = 2.10, (P\ss\P) = 0.08 , (85) 
and therefore, the strange quark content in the proton is rather small, 



y = _ 1 — i-f— = 0.03. (86) 



2(P 


ss 


P) 


(P 


uu 


+ dd\P) 



This value is smaller than the "standard" value y ~ 0.2. However, as we will 
see later, our result depends on the choice of the current u,d-quark mass and 
therefore should not be taken seriously. 
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Figure 16: Dependence of the scalar quark contents in the constituent u-quark on 
the dimension-less coupling constant Cff . The solid line, dashed line and dotted 
line indicate (U\uu\U), (U\dd\U) and (U\ss\U), respectively. 



Assuming the isospin symmetry, the T, nN and T, KN terms are represented 
using the scalar quark contents in the nucleon as follows. 



SttTV = rh (N\uu + dd\N) , 

3 1 - 

Zkn = -(rh + m s )(N\-(uu + dd) + ss\N) 



with rh = {m u + m^)/2. The calculated E^jy an d ^kn terms are shown in Fig. 
17| . The T, nN term is almost independent on G e ^ . At G e p = 0.7, we obtained 
S„-jv — 43.2 MeV which is in good agreement with the value extracted from the 
low-energy 7rN scattering data |j45| : S^tv = (45 ± 10) MeV. As for the T,kn term, 
our result at Gf = 0.7 is 280 MeV. 

In the additive quark assumption of the nucleon states, the following relation 
holds. 

(P\2dd-uu\P) = 3(U\dd\U) . (89) 

Therefore, it vanishes if there exists no flavor mixing in the scalar channel. Our 
result at = 0.7 is (P | 2dd -uu\P) = 0.89. 
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Figure 17: Dependence of the S^tv and E^at terms on the dimension-less coupling 
constant G e § . The solid line indicates T, nN term and the dashed line indicates 
T, KN terms, respectively. 



The pion-nucleon sigma term and the scalar quark contents in the nucleon 
are extensively studied in the three-flavor NJL model in Ref. ||18|| . As pointed 
out there, the scalar quark contents in the constituent quarks depend on the 
current quark masses nonlinearly. In order to study this nonlinearity, we change 
the current u,d-quark mass from 7.5 MeV to 8.5 MeV. Other four parameters of 
the model (m s , Gs Go and A) are determined so as to reproduced the observed 
values of m n , tjik, Tn v and f n . Thus the size of the chiral symmetry breaking is 
fixed though the current u,d-quark mass is changed. The calculated results are 
shown in Table [I]. The scalar quark contents in the constituent u-quark and the 
sigma terms, especially, (U \ss\ U) and T, KN depend on the current u,d-quark 
mass rather strongly. 

We next discuss the validity of the additive quark assumption for the nucleon 
state. Kunihiro and Hatsuda |17|,£|6| studied the effects of the confinement and 
the short-range spin-spin interaction between the constituent quarks in baryons 
in the framework of the nonrelativistic potential model. They have found that 
the effects of the residual interactions between the constituent quarks in the 
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Table 1: The scalar quark contents in the constituent u-quark and the E^jv and 
Y, KN terms in the three different current u,d-quark masses. 



m u4 [MeV] 


(U\uu\ U) 


(U\dd\ U) 


(U | ss | U) 


S.tv [MeV] 


Y, KN [MeV] 


7.5 


1.582 


0.707 


0.210 


51.5 


390.6 


8.0 


1.388 


0.478 


0.091 


44.8 


308.8 


8.5 


1.232 


0.300 


0.009 


39.1 


243.0 



proton increase the uu, dd and ss contents of the proton by about 5%, 24% and 
13%, respectively. The contributions of the kinetic term of the confined quark 
are flavor independent and negative, while those of the short-range spin-spin 
interaction are flavor dependent and positive. The flavor-mixing effect by the 
short-range spin-spin interaction is rather large. For simplicity, we consider the 
(U \dd\U) = (U | ss | U) =0 case. Using the expression of (P \ qiqi \ P) given in 
Ref. |T7|], one obtains 

6b 

(P | 2dd -uu\P) = —(U\uu\U) , (90) 

u 

where b = (176.4MeV) 3 is the strength of the short-range spin-spin interaction. 
Inserting our numerical results at G e § = 0, i.e., M u = 325 MeV and (U \uu\U) = 
1.79, we get (P | 2dd -uu\P) = 1.72. 

Recently, the static properties of the nucleon have been studied in the rela- 



tivistic Faddeev approach using the two-flavor NJL model f47 |. Their results of 
the scalar quark contents in the nucleon are (P | uu \ P) — 1.795 and (P\dd\P) = 
1.095. Since the single quark renormalization factor DM /dm is not included in 
their calculations, one should compare above numbers with those in the additive 
quark model, i.e., (P | uu \ P) = 2 and (P\dd\P) = 1. The residual inter- 
actions between quarks in the nucleon decrease the u-quark content by about 
10% and increase the d-quark content by about 10%. It is a contrast to the 
results in the nonrelativistic quark model. As for the flavor mixing effect, one 
gets (P | 2dd — uu\P) = 0.395([7 \ uu\U) ~ 0.71. Here we have used our re- 
sult: (U | uu | U) — 1.79. In the nonrelativistic quark model approach, the flavor 
mixing effect entirely comes from the short-range spin-spin interaction. On the 
other hand, in the relativistic Faddeev approach, only the scalar diquark state is 
included. Inclusion of the axialvector diquark may be important for the flavor 
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mixing in the scalar quark contents in the nucleon. 

In the above two approaches, the effects of the pion cloud around the quark- 
core of the nucleon is not taken into account. Wakamatsu has studied the 
scalar quark contents in the nucleon using the chiral quark soliton model as 
functions of the constituent u,d-quark mass PpJ. His results follows. 
At M u4 = 350MeV, (P\uu\P) = 1.807, (P\dd\P) = 1.223 and at M u4 = 



450MeV, (P\uu\P) = 1.382, (P\dd\P) = 0.978. As discussed in @, the 



flavor-asymmetry of the sea-quark in the nucleon gives rise to the flavor mix- 
ing phenomena and the calculated results are {P \ 2dd — uu\P) = 0.639 at 
M u4 = 350MeV and (P | 2dd -uu\P) = 0.574 at M u4 = 450MeV. Since the 
chiral quark soliton model do not have the dynamics of the spontaneous breaking 
of chiral symmetry, the effects of the single quark renormalization factor DM/ dm 
is not taken into account. 

In conclusion, as is expected, the flavor mixings in the scalar channel are 
approximately proportional to the strength of the Ua{1) breaking interaction in 
the NJL model. However their magnitudes depend on the choice of the model 
parameters. Furthermore, the scalar quark contents in the proton have many 
origins, and therefore we cannot draw a definite conclusion on the effects of the 
Ua(1) breaking interaction. 

6 Effects of the Ua{1) Anomaly in Baryons 

Since the effects of the ?7a(1) anomaly are rather large in the pseudoscalar 
sector, it is natural to ask if one can see some effects in the baryon sector. It 
was pointed out in p4| that the instanton can play an important role in the 



description of spin-spin forces, particularly for light baryons. The pattern of these 
effects can be very hard to disentangle from one-gluon exchange. The effects of 
the instanton induced interaction in baryon number B = 2 systems were studied 
in | 25[| . It was shown that an attraction between two nucleons is obtained by 
the two-body instanton induced interaction, while the three-body interaction is 
strongly repulsive in the H-dibaryon channel and makes the H-dibaryon almost 
unbound. 

We estimate the effects of the £7a(1) anomaly on the B = 1 and B = 2 
systems by employing the six-quark determinant interaction given in Eq. (|4j) 
whose strength was determined so as to reproduce the observed ?7-meson mass, 
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Table 2: Contribution of the two-body term to octet and decuplet baryons. All 
the entries are in units of MeV. 



wave function 


N T S A 


AS* E* Q 


MIT 
NRQM 


-43.9 -41.2. -41.2 -42.9 
-40.88 -36.6 -36.6 -39.4 


0.12 0.12 
0.07 0.07 



Table 3: Baryon component, SU(3) multiplet, spin, isospin and strangeness of 
the eight channels of two octet baryons. 



channel 


Baryon component 


SU{3) multiplet 


Spin 


Isospin 


Strangeness 


I 


NN 


10* 


1 








II 


NN 


27 





1 





III 


NT 


27 





3/2 


-1 


IV 


NT - NA 


27 





1/2 


-1 


V 


NT, - NA 


10* 


1 


1/2 


-1 


VI 


NT 


10 


1 


3/2 


-1 


VII 


NT - NA 


8 


1 


1/2 


-1 


VIII 


H 


1 








-2 



the 7] — > 77 decay width and the r\ — > 7r°77 decay width, namely, Gp = 0.7. It 
is done by calculating the matrix elements of the the Ua(1) breaking interaction 
hamiltonian with respect to unperturbed states of the MIT bag model and the 
nonrelativistic quark model (NRQM). For B = 2 systems, we only consider the 
(OS) 6 configuration of the six valence quark states. Therefore, the matrix element 
with respect to such a state gives a measure of the contribution of the {7a (1) 
breaking interaction either to the dibaryon or to the short-range part of the 
interaction between two baryons. The determinant interaction induces not only 
three-body but also two-body interactions of valence quarks when the vacuum 
has a nonvanishing quark condensate. The details of the calculation are described 
in 

Table |2| shows the contribution of the two-body term for B = 1. The contri- 
bution to the decuplet baryons vanishes in the SU(3) limit and therefore comes 
only from the SU(3) asymmetry of the quark wave function. The three-body 
term does not contribute to the B = 1 states. Thus the NA mass difference due 
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Table 4: Contributions of the two-body term to the eight channels of two octet 
baryons listed in Table 3. All the entries are in units of MeV. 



wave function 


I 


II 


III 


IV 


V 


VI 


VII 


VIII 


MIT 
NRQM 


-90 
-120 


-85 
-105 


-89 
-102 


-86 
-104 


-94 
-118 


-97 
-117 


-122 
-148 


-163 
-183 



Table 5: Contributions of the three-body term to the H-dibaryon (VIII) and 
strangeness —1 two octet baryon channels (III- VII). All the entries are in units 
of MeV. 



wave function 


III 


IV 




V 


VI 


VII 


VIII 


MIT 


-6 x 10~ 2 


-6 x 10~ 2 


-7 


x 10~ 2 


20.7 


25.1 


40.7 


NRQM 


-5 x 10~ 2 


-5 x 10~ 2 


-5 


x 10~ 2 


28.3 


34.9 


56.1 



to the Ua{1) breaking interaction is about 15% of the observed one. 

We next discuss the case of B = 2. We consider all the possible channels which 
are made of two octet baryons listed in Table ^. Table |] shows the contribution 
of the two-body term. The channel VIII gets the strongest attraction, about 170 
MeV, and the channel VII gets the second strongest attraction. The contributions 
of the three-body term to the H-dibaryon and strangeness —1 channels are given 
in Table [5]. It should be noted that the three-body term has no effect on the NN 
channels, and that the contributions to the channels III, IV and V reflect the 
ST/ (3) breaking in the quark wave function. The contributions of the three-body 
term in channels VI, VII and VIII are remarkable and one will be able to observe 
some effects experimentally. 

We should comment on the difference between the determinant interaction 
used here and the instanton-induced interaction used in ref. ||25|| . The relative 
contributions of the Ua(X) breaking interaction within the baryonic sector or 
within the mesonic sector are similar for the two interactions. However, the ratio 
of those in the baryonic sector to those in the mesonic sector is about |. Namely, 
if one fixes the strength of the interaction so as to give the same mass difference 
of 1] and T]', the effects of the instanton-induced interaction in the baryonic sector 
would be about \ stronger than those of the determinant interaction. After this 
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correction the strength of the present {7a(1) breaking interaction is consistent 
with that used in the calculation of the baryon-baryon interaction in ref. p5 . 



7 Summary and conclusions 

Using an extended three-flavor NJL model that includes the 't Hooft instanton 
induced Ua(1) breaking interaction, we have studied the 77 — > 77 , 77 — > 7/i~/i + , 
77 — > je~e + , and 77 — > 7r°77 decays as well as the 77-meson mass, 77 decay constant 
and flavor SU (3) singlet-octet mixing angle 9 of the 77 meson. Advantages of our 
approach are as follows. (1) Effects of the explicit chiral symmetry breaking by 
the current quark masses and the {7a (1) anomaly can be evaluated consistently 
on the 77-meson decay amplitudes. (2) One can study how the 77-meson properties 
change when the strength of the {7^(1) breaking interaction is changed. 

We have found that the 77-meson mass, the 77 — > 77 , 77 — > 7/i~/i + and 
77 — > 7r°77 decay widths are in good agreement with the experimental values 
when the {7a(1) breaking is strong and the mixing angle 9 is about zero. 

The calculated 77 decay constant is almost independent of the strength of the 
{7a(1) breaking interaction and is close to the pion decay constant. It indicates 
that the 77 meson does not lose the the Nambu-Goldstone boson nature though its 
mass and mixing angle are strongly affected by the {7a (1) breaking interaction. 

Our result of the mixing angle 9 is about zero which is different from the value 
9 ~ —20° obtained in the PCAC + ABJ anomaly approach and ChPT. We have 
discussed the possible origin of this difference in Sect. 4.2. It should be stressed 
here that the 77' — > 77 decay width is used to obtain the mixing angle in the 
PCAC and ChPT approaches. However, since the 77' meson is heavy, it is rather 
questionable to study the 77' — > 77 decay in the PCAC and ChPT approaches. 

Since the 't Hooft instanton induced {7^(1) breaking interaction gives rise to 
the flavor mixing not only in the pseudoscalar qq channels but also in the scalar 
qq channels, we have studied the scalar quark contents in the nucleon, the pion- 
nucleon and the kaon-nucleon sigma terms. The calculated pion-nucleon sigma 
term is almost independent of the strength of the {7a(1) breaking interaction and 
in good agreement with the value extracted from the low-energy ttN scattering 
data. Concerning the flavor mixing effects, we have found that the amount of the 
flavor mixing in the scalar quark contents in the nucleon depends on the strength 
of the spontaneous breaking of chiral symmetry and the residual interactions 
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between the constituent quarks in the nucleon rather strongly. In this sense, the 
scalar quark contents of the nucleon are interesting and important quantities. 
Further studies are necessary. 

We have estimated the effects of the C/a(1) anomaly on the baryon number 
B = 1 and B = 2 systems too. We have found that the NA mass difference 
due to the C/a(1) breaking interaction is about 15% of the observed one and the 
three-body term of the Ua(1) breaking interaction gives 40-50 MeV repulsion in 
the H-dibaryon channel. 

Finally, we should note that the NJL model does not confine quarks. Since the 
Nambu-Goldstone bosons, 7r, K and rj are strongly bound, the NJL can describe 
their properties fairly well. However the r/'-meson state in the NJL model has 
an unphysical decay of rj' — > qq. Therefore, we do not apply our model to the rj' 
meson. In order to study the role of the Ua(1) anomaly on the low-energy QCD 
further, the studies of the //-meson properties are desirable. 
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